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APPROXIMATIONS  TO  DIRECTIVITY  FOR  LINEAR,  PLANAR, 
AND  VOLUMETRIC  APERTURES  AND  ARRAYS 

1.  INTRODUCTION 


BACKGROUND 

Even  after  decades  of  sonar  design,  approximations  to  the  directivity  factor  or  index  of  an 
array  are  sometimes  used  inappropriately.  Many  of  the  approximations  commonly  used  provide 
accurate  directivity  approximations  only  for  the  simplest  of  array  geometries.  As  the  array's  size, 
shape,  weighting,  and  complexity  increases,  there  is  a  renewed  need  for  better  directivity 
approximations. 

Directivity  is  defined  as  the  ratio  of  the  output  signal-to-noise  power  ratio  (SNR)  of  an 
array  to  the  input  SNR  at  an  omnidirectional  element  in  a  spherically  isotropic  noise  field. 
Calculation  of  directivity  requires  integrating  the  magnitude-squared  response  of  the  array  over  all 
angles  of  incidence.  In  spherical  coordinates,  these  arrival  angles  are  denoted  by  an  azimuthal 
angle  6  and  a  polar  angle  (j).  Hence,  calculation  of  the  directivity  requires  a  two-fold  integration 
over  the  angular  space  defined  by  the  azimuthal  and  polar  angles.  For  complex  large-size  arrays 
consisting  of  thousands  of  elements,  directivity  calculations  using  numerical  integration  procedures 
can  be  time  consuming,  even  on  state-of-the-art  computing  systems.  This  report  provides  a 
number  of  approximations  for  estimating  the  directivity  of  linear,  planar,  and  volumetric  arrays, 
which  are  allowed  to  have  arbitrary  real  amplitude  shading  coefficients  and  element  locations.  The 
linear  and  planar  arrays  are  allowed  to  have  directional  element  responses  as  well. 

A  practical  discussion  of  some  of  the  well-known  formulas  for  approximating  directivity  is 
given  by  Bobber',  where  the  directivity  factors  (DFs)  of  three  standard  transducers  zu'e  measured 
and  compared  to  five  theoretical  methods  used  to  calculate  directivity.  Two  of  the  methods  are 
based  on  approximations  by  Molloy^,  which  use  the  simple  formulas  of  DF  =  2  /A  for  a  line 
array  with  an  effective  aperture  length  of  L^,  and  DF  =  AtcAJ}}  for  a  baffled  planar  array  with  an 
effective  acoustic  aperture  area  of  A^.  For  both  cases,  the  effective  aperture  includes  the  physical 
array  size  plus  half  the  array's  interelement  spacing  at  each  end.  In  the  free  field,  the  planar  array’s 
directivity  factor  would  be  reduced  by  one-half. 


These  earlier  approximations  assumed  an  ideal  model  for  the  array  response,  where  the 
mainlobe  was  relatively  narrow  (the  length  of  the  aperture  exceeding  several  wavelengths)  and  the 
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peak  sidelobe  levels  were  much  less  than  the  mainlobe  level.  Furthermore,  it  was  assumed  that  flat 
amplitude  shading  weights  were  applied  to  the  individual  array  elements.  Bobber  does  provide  a 
short  discussion  on  shading-weight  correction  factors  that  can  be  used,  in  some  cases,  to  adjust  the 
directivity  factors.  However,  with  current  array  configurations  and  signal  processing,  these  earlier 
approximations  are,  in  general,  too  limited. 

GENERAL  MODEL  FOR  ARRAY  RESPONSE 

For  an  array  with  arbitrary  geometry  and  element  locations,  the  beamformed  array  single¬ 
frequency  amplitude  response  may  be  written  as 

A(e,0)  =  FJ,  (1) 

where  A{6,(p)  is  the  array's  angular  response  to  an  arrival  from  azimuthal  angle  G  and  polar  angle 
0;  N  is  the  total  number  of  array  elements;  and  g„(0,0)  are,  respectively,  the  amplitude 

shading  and  element  amplitude  angular  sensitivity  of  the  n-th  array  element;  k  is  the  wavevector 
corresponding  to  an  arriving  acoustic  planewave  from  (6,<p);  and  is  the  wavevector  to  which 

the  array  is  steered.  The  corresponding  steering  angles  are  0^,0^.  The  «-th  array  element  position 
is  defined  by  the  coordinate  vector  F„  =(x„,y^,z„). 

The  linear,  planar,  and  volumetric  array  coordinates  are  given  in  figure  1 ,  along  with  the 
definition  of  spherical  angles  0,0. 
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z, 


(a)  Linear  (b)  Planar 


z 


(c)  Volumetric 


Figure  1.  Coordinate  Systems  for  the  Linear,  Planar,  and 

Volumetric  Arrays 


3/(4  blank) 


2.  APPROXIMATIONS  TO  DIRECTIVITY 


The  directivity  index  (DI)  is  defined  as  a  decibel  measure  of  the  improvement  in  the  SNR 
that  a  beamformed  array  provides  in  an  ideal  isotropic  noise  field  with  a  perfectly  correlated  signal, 
relative  to  an  omnidirectional  array  element  in  the  free  field,  i.e., 

Dl  =  10  log  DF  ^  10  Log  ,  (2) 


where  DF  is  the  directivity  factor. 

With  the  aid  of  equations  (1)  and  (2),  the  components  that  comprise  the  calculation  of  the 
directivity  index  can  be  written  as 


5^  .  =  5, 

jj—omm 

(3) 

s  =sA(e,(t>) 

array  ^  s 

2 

(4) 

In  n 

N,-^.  =  1  jA!  sin(«#<ie, 

0  0 

(5) 

27LK 

(6) 

0  0 


where  S  is  defined  as  the  received  signal  power  level  and  N-  is  defined  as  the  isotropic  noise  power 
level. 


Noting  that  equation  (5)  is  simply  4;riV,  and  canceling  the  signal  and  isotropic  noise  power 
levels,  the  DF  in  equation  (2)  becomes 


Ikiz 

j  J \A{6,  0)|^  s,m{^)d^d6 


(7) 
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where  |yl(0j,0j|'is  the  array  power  response  at  steering  angle 

Equations  (6)  and  (7)  indicate  an  integral  with  limits  that  account  for  all  acoustic  arrival 
angles.  Hence,  the  directivity  expression  (equation  (7))  assumes  an  unbaffled,  or  free-field,  array 
configuration.  In  a  baffled  configuration,  would  be  reduced,  resulting  in  an  increased 
directivity. 

The  array  amplitude  response  A{d,(p),  given  by  equation  (1),  is  for  the  most  general  three- 
dimensional  array  with  arbitrary  element  locations,  weights,  and  element  responses.  For  this  most 
general  case,  the  sum  of  N  terms  can  be  evaluated  only  by  calculating  every  individual  complex 
term  and  summing  up  those  N  quantities.  This  is  a  particularly  time-consuming  task,  especially 
when  the  number  of  elements  N  is  large.  When  coupled  with  the  fact  that  the  DF  requires 
computation  of  a  double  integral,  namely  equation  (7),  which  must  be  repeated  for  each  different 
frequency  and/or  steering  angle,  the  computational  burden  often  becomes  excessive.  Furthermore, 
as  the  number  N  of  elements  increases,  the  array  power  response  | A(0, 0)|  becomes  even  sharper 
in  angles  6,<p,  thereby  requiring  still  finer  evaluation  of  the  integrand  in  equation  (7)  in  order  to 
retain  accuracy  in  the  final  DF  calculation.  Hence,  the  need  for  accurate  directivity  approximations 
is  obvious. 

For  the  most  general  array  configurations  and  processing,  accurate  directivity  predictions 
can  only  be  obtained  by  numerically  integrating'  the  double  integral  given  in  equation  (7).  Also, 
for  certain  arrays,  the  maximum  array  power  response  ^A{9^ ,  )|  must  be  calculated  and  cannot 

be  estimated.  The  maximum  array  response  may  not  occur  precisely  at  the  steering  angle  (6^  (p^) 
because  of  array  curvature,  unequal  array  element  spacing,  variable  element  sensitivities,  and 
shading  weights.  This  report,  however,  considers  a  broad  set  of  array  geometries  and 
configurations  which  may  be  approximated,  thereby  avoiding  the  need  for  numerical  integration. 
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DIRECTIVITY  FACTORS  FOR  UNIFORM  LINEAR  ARRAYS 


Approximating  directivity  plainly  requires  approximating  the  double  integral  in  equation 
(7);  the  numerator  typically  can  be  determined  directly  from  equation  (1).  In  fact,  for  a  line  array 
with  constant  interelement  spacing,  the  double  integral  in  equation  (7)  may  be  evaluated  exactly  at 
the  frequency  for  which  spacing  dy  =  A/2,  regardless  of  the  steering  direction.  This  leads  to 

Z)/  =  101og|(^^^^^w„)  independent  of 

and  for  the  special  case  of  flat  weighting,  to  the  result 
DI  =  101og{A/}  =  101og{24/A} 

for  d,  -  A/2,  as  mentioned  in  section  1 .  Although  this  is  an  exact  result  for  the  uniform  linear 
array,  it  is  limited.  For  example,  figure  2  illustrates  the  exact  DF,  obtained  by  numerical 
integration  and  normalized  by  the  number  of  array  elements,  for  an  equispaced  uniformly-weighted 
line  array  with  nondirectional  elements  (uniform).  For  the  frequency  at  which  /  =  /„,  where 
/c(=  cl^d^  )  is  the  design  frequency,  the  DF  is  independent  of  array  steering  and  equals  N.  This 
is  valid  only  at  the  design  frequency  /^;  as  the  line  array's  operating  frequency  changes,  the  DF 
becomes  a  function  of  frequency  and  steering  angle.  When  /</,,,  a  flare  in  the  DF  occurs  near 
endfire,  cos(^^)  ^  ±1,  and  results  in  a  doubling  of  the  DF  precisely  at  endfire  steering.  This 
doubling  results  from  the  elimination  of  all  grating  lobes  within  the  array  response  region. 
Equivalently,  at  these  frequencies,  the  array  response  may  be  considered  oversampled  or 
overpopulated,  J,/A  <  l/2.  Conversely,  a  droop  in  the  DF  occurs  near  endfire  steering 
when  /  >  because  of  the  presence  of  grating  lobes  within  the  array  response  region.  As  the 

operating  frequency  increases  above  the  design  frequency,  the  droop  occurs  further  from  endfire; 
in  fact,  the  breakpoint  is  approximately  where  cos((^^.)  =  2/^/ /  -1-l/V. 

Derivations  for  the  DF  of  more  general  linear  arrays  are  presented  in  detail  in  the  following 
pages.  Subsequent  approximations  (for  planar  and  volumetric  arrays)  follow  the  methodology 
outlined  here  and  are  less  elaborate  in  their  presentations. 
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2.0 


/V  =  40 


0.0  0.25  0.5  0.75  1.0 

Us  =  COS  ^s 


Figure  2.  Normalized  Directivity  Factors  for  a  Uniformly  Weighted 

Equispaced  Line  Array 


DISCRETE  LINEAR  ARRAYS 

The  line  array  is  oriented  along  the  z-axis,  as  shown  in  figure  1,  and  each  array  element  is 
allowed  to  be  directional.  For  example,  the  elements  can  have  a  common  cosine  angular 
dependence.  The  maximum  angular  response,  or  boresight,  of  these  elements  is  chosen  to  be 
collinear  to  the  y-axis.  Therefore,  the  amplitude  angular  sensitivity  becomes 

0)  =  |cos(0)[  =  |sin(0)sin(0)|  ^  where  ©  =  (0, 0)  is  the  conical  angle  between  the  element's 
boresight  and  the  planewave  arrival  angle,  and  v  is  an  arbitrary  power  to  which  the  cosine  is 
raised.  For  omnidirectional  elements,  v=  0.  Figure  3  illustrates  some  typical  element 
directionality  patterns. 
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V  =  0.0 


Figure  3.  Element  Angular  Sensitivity  |cos(0)|  for  Powers  v  =  0,  0.5,  1 

For  the  line  array  geometry,  each  element  has  the  same  angular  dependence;  that  is, 
is  independent  of  n.  Hence,  from  equation  (1),  the  double  integral  in  equation  (7) 

becomes 


2k  n 


0  0 


^  =  J  1  Xli  w«|sin(0)sin(0)[  exp{/fc:„(cos(0)  -  cos((^J)} 


sm{(j))d([)dd . 


(8) 


Setting  u  =  cos(^)  and  =  cos(^^)  yields 


2k  \  ^ 

/  =  J {sin^(0)j'' J {l  -  exp{zfc„(w  -  mJ} 


-1 


du. 


(9) 


For  integer  v,  the  integral  over  9  is  straightforward;  in  fact,  noninteger  vcan  be  evaluated 
using  gamma  functions/  Thus,  there  follows 
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(10) 


1 

I  =  /2(v)  j (1  -  u-y  -  wj} 


-1 


du. 


where 


h{v)  = 


2<^^"'>r-(v  +  l/2) 
r(2v+i) 


=  2-\Ik 


r(v4-i/2) 

r(v+i)  • 


(11) 


The  simplification  is  obtained  by  the  use  of  the  duplication  formula^  for  gamma  functions.  If  it  is 
assumed  that  the  elements  have  a  cosine  angular  dependence,  then  v  =  1  and  h{l)  =  it,  whereas 
for  omnidirectional  elements,  v=  0  and  h(0)  =  2k. 


The  focus  here  is  on  the  behavior  of  the  sum  in  equation  (10)  for  large  N.  Using 
Lagrange’s  trigonometric  identity,  it  is  easily  shown  that  for  uniform  weighting  and  equidistant 
element  spacing  d,  ,  the  absolute  value  of  the  series  sums  to 


exp^ikndyu  -  wjj  =  |sin(A^M_.(w 


)l2)/sm{kdyu  -  u^)/2) 


(12) 


The  square  of  this  sum,  normalized  by  its  maximum  value  N,  is  plotted  in  figure  4  along  with  the 
slowly  varying  m -function  1  -  u'  over  the  integration  limits  given  in  equation  (10),  for  the  case 
where  c/.  =  A/2  and  =  0. 

It  is  apparent  that  the  major  contribution  to  the  integral  (10)  occurs  at  u  =  u^  and,  provided 
that  =  cos(0J  is  not  close  to  ±  1  and  that  <  A/2 ,  the  contribution  of  the  integrand  function 
1  -  u'  may  be  approximated  at  u  =  u^.  The  integration  in  equation  (10)  then  simplifies  (for  real 

weights)  to 


1 

I  =  h{v){\ -  u^y  j  exp{/a;r(n  -  m){u  -  u^^du,  (13) 

-1 


where  or  is  a  number  proportional  to  half-wavelength  spacing,  i.e.,  d,  =  a  A/2.  By  interchanging 
the  order  of  the  u -integration  with  the  summations,  the  integral  can  be  evaluated  exactly  as 
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The  Kronecker  delta  is  denoted  (5„„.  For  half-wavelength  element  spacing,  a  =  l  and  equation 
(14)  simplifies  to 

/s2/i(v)(sin((^J)'''Xlvv„l  (15) 

The  maximum  array  response,  in  the  direction  the  line  array  is  steered,  can  be  evaluated 
directly  from  equation  (1)  as 

=sin^''(6>,)sin^’'((/»,)|^^^jW„  .  (16) 


u 

Figure  4.  Comparison  of  the  Slowly  Varying  Function  1-  u^  to  the  Function 

K  K  ^  ^ 

%m{—Nu)  Nsm{—u)  for  =  2,  5,  20 

2  !  \  2  ) 
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Therefore,  substitution  of  equations  (14)  and  (16)  into  equation  (7)  yields  the  approximation  for 
the  DF  of  a  linear  array  of  equispaced  directional  elements  having  real,  but  arbitrary,  shading 
weights,  namely. 


m 


2 

hi  M 

>  > 

|<5„„+(l-5„„)cos(a;r(n- 

m)cos(</i,)  -i- - ^ 

anyn-m) 

(17) 


The  above  expression  could  be  written  in  terms  of  frequency  using  the  substitution  a  =  2d,J jc, 
with  /  the  operating  frequency  in  Hertz  and  c  the  sound  speed,  provided  that  /  <  c/(2J,)  =  f^. 


For  cosine  element  sensitivity  (v  =  1)  and  half-wavelength  element  spacing,  the  DF  in 
equation  (17)  simplifies  to 


DF  =  - 


2sin"(6Jj 


N 


s2sin“(e,)W„ 


(18) 


where  is  the  effective  number  of  weights.  For  the  array  element's  maximum  response,  at 
6^  =±  njl  (in  the  y-axis  direction),  the  cosine  element  sensitivity  provides  a  doubling  of  the  DF 
compared  to  that  of  an  omnidirectional  element.  Thus,  cosine  element  directivity  can  provide,  at 
broadside,  a  3-dB  increase  in  the  line  array's  DI. 

Approximations  for  discrete  linear  arrays  with  nonperiodic  element  spacing  (provided  no 
interstitial  array  spacing  is  greater  than  half  an  acoustic  wavelength)  are  derived  in  the  following 
section  as  a  special  case  of  a  weighted  continuous  linear  aperture. 
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CONTINUOUS  LINEAR  APERTURES 


For  a  continuously  weighted  linear  aperture,  the  amplitude  response  for  steering  angle 
can  be  written  as 

A{d,(t>)=  jw(z)g(z,0,0)exp{fc(cos(0)-cos(0,))}ifz,  (19) 


where  the  real  weighting  w{z)  automatically  terminates  the  integral  at  finite  limits.  This  expression 
corresponds  to  the  discrete  array  response  shown  in  a  general  form  by  equation  (1).  Assuming,  as 
in  the  discrete  case,  that  this  aperture  has  an  amplitude  angular  sensitivity  g(6,  (p)  (independent  of 
z),  which  is  proportional  to  a  power  of  the  cosine  of  the  angle  between  the  aperture's  maximum 
response  axis  and  the  planewave  arrival  angle,  then 


A{d,<p)  =  |sin(0)sin(^)|  |  w(z)exp|zAz(cos(^)  -  cos(<pJ)^dz. 


(20) 


For  uniform  angular  response,  the  power  is  set  to  zero,  v  =  0. 

Substituting  equation  (20)  into  the  double  integral  in  equation  (7)  and  letting  u  =  cos{(p)  as 
before,  yields,  after  some  simplification, 


/  =  /i(v)J  j  w{z)txp{ikz(u - u^}dz 


du , 


(21) 


where  h(y)  is  defined  by  equation  (11).  Three  steps  are  necessary  to  approximate  this  integral. 
First,  as  done  in  the  discrete  linear  case,  it  is  argued  that  (1  -  u^Y  varies  slowly  with  wand  that  its 
contribution  is  well  approximated  by  its  value  at  u  =  u^.  Thus,  for  |wj|  =  |cos(0j|  <  1, 


/s/i(v)(l 


|2 


w{z)e'^'dz 


du, 


(22) 


where  T  =  k(u  -uj.  The  second  step  for  the  integral  approximation  considers  the  Fourier 
transform  of  the  weighting  function  w(z),  which,  by  definition,  is 
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oo 


(23) 


W{t)  =  jw(z)e‘^dz. 


If  it  is  presumed  that  the  weighting  >v(z)  is  even  in  z,  then  a  Gaussian  approximation  for  W(t), 
which  can  be  applied  for  all  T,  is* 


lV(T)  =  fV(0)exp 


T' 

—  m 

V  2  , 


(24) 


where  the  mean-square  length  of  weighting  w(z)  is  defined  as 


f  z^w(z)dz 

^ - .  (24a) 

I  w{z)dz 


This  choice  matches  the  moments  of  >v(z)  through  second  order.  Therefore,  recognizing  that  the 
u  -integration  limits  may  now  be  extended  to  ±  <» ,  equation  (22)  becomes 


/  =  /2(v)|sin"(0^)|  IT^(O)  Jexp{-A:^(w  -  u^fm^}du. 


(25) 


Finally,  the  integral  in  equation  (25)  may  be  evaluated  exactly,  yielding 


/  =  A(v){sin^(0,)}''w-(O) 


(26) 


Thus,  the  approximation  for  the  DF  of  a  continuous  line  aperture,  which  is  allowed  to  have 
nonuniform  even  weighting  and  a  directional  angular  response,  is 


h(v)  A  A 


V^r(v  +  l)sin-'^(e,)' 
r(v  +  l/2) 


(27) 


*An  alternative  approach  for  a  linear  aperture  is  given  in  apptendix  A. 
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The  quantity  m  = 


|j_  z^w(z)dz 


w{z)dz 


1/2 

in  equation  (27)  may  be  interpreted  as  the  root- 


mean-square  (RMS)  length  of  weighting  w(z). 


For  maximum  directivity,  the  best  value  for  steering  angle  9^  is  9^  =  ±;r/2  (along  the  y- 
axis).  It  is  illustrative  to  consider  the  v-dependence  of  the  bracketed  term  (shading  factor)  in 
equation  (27)  on  directivity.  Recall  that  as  v  increases,  the  aperture’s  angular  response  becomes 
narrower.  Figure  5  plots  the  bracketed  term  in  equation  (27)  for  various  fixed  values  of  steering 
angle  9^.  Notice  that  for  the  steering  angle  9^  =  kI1,  the  function  starts  at  unity  for  v  =  0  and 

increases  monotonically  with  v.  Given  this  choice  of  steering  angle,  it  is  seen  that  the  DF 
increases  with  element  directivity.  However,  any  other  choice  of  steering  angle  leads  to  an 
eventual  decay  of  the  DF  with  v.  This  implies  that  in  general,  highly  directional  elements  will 
result  in  a  gain  in  directivity  along  the  element's  boresight,  but  with  the  cost  of  a  significant 
decrease  in  directivity  for  other  steering  angles. 


V 


Figure  5.  Shading  Factor  for  Line  Aperture  with  jcos(©|  Element  Directionality 

It  should  be  noted  that  the  DF  in  equation  (27)  is  independent  of  polar  steering  angle 
and,  hence,  of  sin(0J.  More  generally,  it  may  be  shown  that  if  the  element  amplitude  response 
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g(z,6,(l>)  in  equation  (19)  varies  slowly  with  angle  and  is  independent  of  position  z,  then  equation 
(27)  is  replaced  by 


^  \-J g\9,<p,)de 


(27a) 


The  latter  factor  is  obviously  unity  for  a  flat  angular  response  g(6,<j>).  If  the  angular  response  is 
separable,  that  is,  g(6,(l))  =  g^id)g2{(l>),  the  approximate  DF  simplifies  to 


(27b) 


which  is  independent  of  the  polar  angular  dependence  g2(<t>)-  This  independence  of  the  DF  on  the 
element  directivity  will  also  appear  in  the  derivations  for  the  planar  and  volumetric  apertures.  The 
special  result  in  equation  (27)  corresponds  to  ^,(6)  =|  sin(e)  T  and  any  slowly  varying  gji^)- 
For  a  uniform  angular  aperture  response,  v  =  0,  equation  (27)  simplifies  to 

DF  =  a4k  —  .  (28) 

A 


From  equation  (27)  or  equation  (28),  it  is  straightforward  to  obtain  an  approximation  for 
the  directivity  of  &  discrete  line  array  having  arbitrary  element  locations  {z„},  provided  that  no 

interstitial  element  spacing  is  greater  than  a  half-wavelength  and  the  origin  of  coordinates  is  located 
so  that  21^-1  example,  a  discrete  line  array  of  N  elements  has  an  impulsive 

weighting  function 

w{z)  =  Y^w„5{z-z„).  (29) 


Substituting  this  weighting  function  into  equations  (24a)  and  (28)  gives 


DF 


_  U 


1/2 


■  {i-r 


(30) 
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DIRECTIVITY  FACTORS  FOR  PLANAR  APERTURES 


The  approximation  ATtAg/X^  for  the  DF  of  an  equispaced  baffled  planar  array,  where  Ae  is 
the  effective  array  area,  is  only  applicable  for  an  array  having  ideal  half-omnidirectional  elements, 
flat  weighting,  and  broadside  steering.  These  restrictions  can  be  eliminated  by  using  the  Gaussian 
approximation  technique  described  above  for  the  continuous  linear  aperture. 


The  amplitude  response  for  a  continuous  planar  aperture  with  arbitrary  real  weighting 
(which  need  not  be  separable)  and  lying  in  the  xy-plane  is 

=  |cos(^)|^  [  J  w{x,y)exp{iksm(<l))[xcos(6)  + ysm(6)] 

(31) 

-  iksin(^^)[xcos{6J  +  ysin(d^)]}cixdy, 


where  it  is  presumed  that  the  aperture's  amplitude  directionality  is  proportional  to  a  cosine  to  a 
power,  with  the  maximum  response  or  boresight  aligned  along  the  z-axis,  (f>=0,  namely, 
broadside  to  the  aperture.  The  geometry  is  illustrated  in  figure  lb.  A  uniform  or  omnidirectional 
aperture  angular  response  corresponds  to  v  =  0 .  Again,  the  weighting  function  H’(x,y)  terminates 
the  integral  and  defines  an  effective  finite  aperture.  Equation  (31)  may  be  written  more  compactly 
by  introducing  the  multidimensional  Fourier  transform 

W(u,v)  —  j  J  M>{x,y)Qxy>{iiix +  ivy}dxdy .  (32) 


Now  let  u  =  sin(0)cos(0)  and  v  =  sin(^)sin(0) .  Then,  the  planar  aperture  amplitude 
response  becomes 

A(8,<t>)  -  |cos((^)['^PF(^w  -  ku^,kv  -  kvj_ 

Evaluation  of  the  double  integral  in  equation  (7)  requires  a  one-to-one  mapping  of  the  0^plane  into 
the  Mv-plane.  The  Jacobian  of  the  transformation  is 

=  cos(<^)  sin((^>) .  (33) 

<7(0,0) 
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To  prevent  a  two-to-one  mapping,  the  upper  limit  on  the  0 -integral  in  equation  (7)  must  be 
decreased  to  <p2  =  njl .  This  change  in  the  limit  is  possible  because  the  planar  aperture  response  is 
symmetric  about  the  :xy-plane  (see  equation  (31)).  This  observation  allows  for  the  change  in  the 
upper  limit,  and  the  0 -integral  is  simply  doubled;  if  the  integral  over  (0,:;r/2)  is  not  multiplied  by 
two,  the  resultant  DF  will  be  for  an  ideally  baffled  planar  array.  With  this  change,  the  integration 
region  in  the  0^-plane  maps  to  the  interior  of  a  unit  circle  Q  in  the  wv-plane.  Thus,  from 
equations  (31)  to  (33),  the  denominator  of  equation  (7)  becomes 

1  =  ijj (l-u^  -  v")''"'^'  |^(A:w  -  ku^,kv  -  A:v dudv.  (34) 

c, 

The  peak  of  the  integrand  in  the  Mv-plane  occurs  at  namely  at  radius 

+  v]  =  If  sin(^^)  <  1  ,  that  is,  (p^  <  ;r/2  or  steered  to  the  upper  hemisphere,  then 

cos(0J  >  0  and  the  integral  may  be  approximated  in  the  manner  described  earlier,  namely, 

^  ^  2(co£(^jy^j~  ^°‘J^I^2i,v)\dudv.  (35) 

k 


The  Gaussian  approximation  for  the  double  Fourier  transform  of  the  (even)  weighting 
function  w(x,y)  (as  in  equation  (24)  for  the  continuous  line  aperture)  is* 


f 

W(u,v)  =  W[0,0)exp 

V 


(36) 


where  mean  square  lengths 


7 


J  J  x'w{x,y)dxdy 

Aoo  m 

J  I  w(x,y)dxdy 


w(x,y)dxdy 


(36a) 


Substituting  the  Gaussian  approximation  into  equation  (35)  and  performing  the  uv-integration 
yields 


'An  alternative  approach  for  a  planar  aperture  is  given  in  appendix  A. 
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(37) 


^  _icos-i<l>,)y-^/-W\0,0)X^ 
iTOnm^ 

for  the  denominator  of  equation  (7),  where  it  has  been  assumed  that  »  A  and  m^.  »  A . 

Hence,  the  DF  for  an  unbaffled  planar  aperture  is 

DF  =  8;r^cos(0J-^  (38) 

A 

for  which  is  independent  of  power  v.  The  quantity  may  be  interpreted  as  the 

RMS  aperture  area  in  the  j:y-plane. 

More  generally,  if  the  factor  |  cos(0)  |''  in  equation  (31)  is  replaced  with  a  slow  angular 
dependence  g{d,  </>),  as  in  equation  (19),  it  may  be  shown  that  exactly  the  same  result  for  the  DF 

(equation  (38))  is  obtained.  The  approximate  planar  DF  is  independent  of  a  general  slowly  varying 
angular  response  g(d,<p).  The  cos(0^)  dependency  in  equation  (38)  results  from  the 
foreshortening  of  the  planar  aperture  in  the  xy-plane  when  seen  from  the  polar  angle  0^ ;  it  has 

nothing  to  do  with  the  aperture's  angular  sensitivity.  The  reason  for  the  independence  on  angular 
sensitivity,  g(6,(p),  or  on  parameter  v,  for  a  planar  aperture  is  due  to  its  ability  to  resolve  in  both 
6  and  0.  The  narrow  resolution  of  the  planar  aperture  response  ^4(0,0)  quickly  overcomes  that 
of  the  relatively  broad  and  slowly  varying  angular  response  g{9, 0) ;  this  assumes,  of  course,  that 
V  is  of  moderate  order  for  dependency  |cos(0)|  . 

As  an  example  of  approximating  the  directivity,  consider  the  simple  case  of  a 
nondirectional,  uniformly  weighted  rectangular  aperture  steered  to  broadside  ( 0^  =  0).  Letting 
and  be  the  physical  lengths  of  the  aperture  in  the  x-  and  y-directions,  respectively,  then 
and  m^.  =  L^.J (2^3).  Thus,  the  DF  becomes 

DF  =  2k-^  =  -2k-^.  (39) 

3A'  3  A' 

The  well-known  approximation  for  an  unbaffled  planar  array  steered  to  broadside  and  having 
uniform  weights  is  DF  s  2tiAJ 2} ,  where  =  L^L^.  (see  appendix  A).  Therefore,  the  above 
Gaussian  approximation  is  5  percent  larger,  giving  roughly  0.2-dB  greater  DI. 
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Extending  the  Gaussian  approximation  for  directivity  (equation  (38)),  to  discrete 
equi  spaced  planar  arrays  presents  no  complications,  provided  that  the  array  spacings  are  not  greater 
than  half  an  acoustic  wavelength  <  Xjl .  For  instance,  the  weighting  function  for  an 

equispaced  rectangular  planar  array  may  be  expressed  as 


N  M 


=  x-d^ 


^=1  m  =  l 


A^  +  1 


n  ■ 


q  y-<. 


jj 


M+i 


m  ■ 


(40) 


where  need  not  be  separable  in  n  and  m.  From  equations  (36a)  and  (38),  the  directivity 
factor  for  the  discrete  equispaced  unbaffled  planar  array  follows  immediately  as 


N  M 


w. 


iv+r 


2-11/2 


N  M 

IS 

I  n  =  l  m=l 


M  +  1 


m  ■ 


N  M 

IS 

n=l  m=l 


(41) 


for  (p^  <kI1. 

For  flat  weights,  expression  (41)  simplifies  to  the  compact  form 

K  dd^,  i  ,  \l/2/  „  n1/2 

Z)F  =  -2;rcos((^J-^  A/--1)  (M'-i)  . 

If  d^  [N'  -  is  interpreted  as  the  RMS  length  of  the  array  in  the  x  direction  (with  a 

similar  interpretation  in  the  y  direction),  this  expression  for  DF  simplifies  to  , 

which  is  identical  to  the  DF  given  for  the  unbaffled  planar  aperture  in  equation  (38).  The  quantity 
behaves  as  A/ - 1/(2^)  for  large  TV,  leading  to  ~  d^  N 1 . 

More  generally,  for  arbitrary  element  locations,  the  weighting  function  is  given  by 

T 

w{x,y)  =  '^Wj5(x  -  x,)5(y  -  y,),  where  T  is  the  total  number  of  elements  and  locations 

r=l 

are  arbitrary,  except  that  no  interstitial  array  spacings  can  exceed  a  half-wavelength. 


20 


DIRECTIVITY  FACTORS  FOR  VOLUMETRIC  APERTURES 


Exact  expressions  for  the  directivity  of  a  solid  spherical  volumetric  aperture  and  a  spherical 
shell  are  given  by  Anderson/  For  values  of  Z)/A  » 1,  where  D  is  the  diameter  of  the  sphere,  the 
expressions  for  the  DFs  yield  the  approximations 


A, 


and 


DF,,,,  =  U.6\ 


fD\ 

U 


(42) 


(43) 


These  expressions  are  for  continuously  distributed  elements  that  are  uniformly  weighted  and 
nondirectional.  The  exact  DF  for  a  thick  shell  of  inner  and  outer  radii  of  r,  and  r2,  respectively, 
follows  readily  from  the  procedure  in  reference  (6),  although  not  done  there.  Here,  the  authors’ 
approximation  procedure  (to  be  derived)  yields  the  form 


1  +  rjV/  +  r,  r,'’  +  r/ 

5A^  r2  +  r,  r,  +  r/ 


(44) 


For  a  solid  spherical  aperture,  r,  =  0,  and  for  a  shell,  r,  =  Tj.  Thus,  the  authors’  Gaussian 
approximations  are 


and 


An- 
5  A' 


4;r^ 

3  A' 


13.16^. 

A' 


(45) 

(46) 


These  approximations  compare  reasonably  well  with  equations  (42)  and  (43),  respectively. 

The  approximate  DF  derived  below  allows  for  arbitrary  volumetric  shapes,  element 
positions,  and  shading  weights;  however,  the  elements  must  remain  nondirectional.  For  certain 
volumetric  geometries,  such  as  a  cylinder  or  sphere,  it  may  be  possible  to  allow  element 
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directionality,  although  the  derivations  are  likely  to  be  extremely  difficult.  Furthermore,  a  baffle 
may  be  introduced  in  the  volumetric  aperture  so  that  steerings  and  planewave  arrival  directions  are 
restricted  to  the  upper  hemisphere  (see  figure  Ic);  the  volumetric  aperture  does  not,  in  general, 
have  a  symmetric  response  about  the  xy-plane. 

The  derivation  of  the  approximate  DF  here  follows  closely  the  Gaussian  approximation  for 
the  continuous  planar  aperture.  The  amplitude  response  of  a  general  continuous  volumetric 
aperture  is 

A(0,(p)=\  J  I  w {x,y,z) exp {ik[x  cos{6)  sm{(p)  +  ysm{d)sm{4>)  + zcos((())] 

(47) 

-  zA:[xcos(0Jsin(0J  +  >’sin(0Jsin(</>J  +  zcos((pJ]}dxdydz, 

where  w(x,y,z)  is  a  general  three-dimensional  continuous  real  weighting  function.  Notice  that 
the  angular  response  of  the  aperture  itself  is  assumed  omnidirectional.  Defining  the  triple  Fourier 
transform 

W(u,v,w)  =  j  J"  J  w(x,y,z)exp{iux  +  ivy  +  iwz}dxdydz  (48) 


for  all  u,  V ,  and  w,  allows  the  amplitude  response  to  be  recast  in  the  form 

A(6,(p)  -  W (ksm(4>)cos(9)-  ku^,ksin((l>)sm(d)-  kv^,kcos((p)  —  kw^)  (49) 

for  0  <  0  <  27t,  0<  (p  <7t,  where  -  sin(0Jcos(0^),  =  sm(^^)sin(0^),  and  =  cos(0J. 
Note  that  the  amplitude  response  in  the  steered  direction  (0^,  0^)  is  simply 

/i(0„0,)  =  jr (0,0,0)  =  w(x,y,z)dxdydz.  (50) 

Substituting  equation  (49)  into  the  double  integral  in  equation  (7)  and  employing  the  same 
variable  transform  and  arguments  used  to  develop  equations  (32)  to  (34)  gives,  for  the 
denominator  of  equation  (7), 


22 


(51) 


IJ 


dudv 


2  2 

U  -V 


wiku-  ku^,kv  —  kv^,k4\  —  u~  —  v~ 


Note  that  equation  (51),  unlike  equation  (34),  is  not  multiplied  by  2  and  thus  corresponds  to  a 
baffled  volumetric  aperture,  which  responds  only  in  the  upper  hemisphere,  (/»  <  ;r/2. 

A  new  obstacle  for  the  volumetric  aperture  is  the  quantity  ,  which  appears 

in  W .  This  quantity  would  become  imaginary  for  u~  >\  \  hence,  a  small-order  expansion 
about  and  is  introduced.  That  is,  let  u  =  +  £  and  v  =  -f  5-  Equation  (51)  then  becomes 


_ d^ _ 

7i-(w,  +  £)‘-(^  +  ^)- 


\W{ke,k5,kQ{e,5)f, 


(52) 


where  Q{£,  5)  =  '\/l  -  (w^  4-  ef  -  (v^  -f-  df  -  ^  1  -  «/  -  v/  and  C,  is  a  unit  circle  centered  at 
{£,5)  =  {-u^,-vf) .  The  term  Q{£,5)  can  be  linearly  approximated  by  invoking  a  binomial 
expansion  and  neglecting  terms  of  second  order  in  £  and  5 .  Thus,  Q{£,6)  =  —(u^£  +  v^5)/w^ 
for  small  £,5. 


For  symmetric  and  real  weighting  function  w(x,y,z),  the  Gaussian  approximation  of 
W  {u,v,w)  for  all  u,  v,  and  w  (as  in  equation  (36)  for  a  planar  aperture)  is  written  as 


W(u,v,w)  =  Wi0,0,0)exp 


w 


.2  ^ 


- - m'. - m 

y  2  ^  2  ^  2  J 


(53) 


where 


ml  = 


and 


ml  = 


[  J  J  x^w{x,y,z)dxdydz  ^  111  y^^{x,y,z)dxdydz 

_  ^ - 

moo  y  ^oo  moo  moo  ? 

w(x,y,z)dxdydz  J  J  J  '^(x,y,z)dxdydz 

I  J  J  z^w(x,y,z)dxdydz 
_ 

^w(x,y,z)dxdydz 


(53a) 
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Using  this  Gaussian  approximation,  it  follows  that 


\W  (ke,  kS,  kO{£,  (5))|‘  0, 0)exp 


—k'S^ml  -  k'5~m\  -  k 


2  {eu^  +  5v;f 


m: 


(54) 


for  all  £,5.  Substituting  equation  (54)  into  equation  (52),  and  recognizing  that  the  limits  of 
integration  may  now  be  extended  to  ±°° ,  yields 


/-  — r  rr'(0,0,0)exp 


a  2 

- £ 

2 


—5'  +7^(5 
2 


d£d5. 


(55) 


where  a  =  2k^{ml+mlul  ^  =  2k'{m'^+m:v~Jw]^,zn6.  J  - The 
integrals  in  equation  (55)  may  be  evaluated  exactly'*,  provided  that  a  >  0 ,  jS  >  0,  and  aP>  ; 
thus,  the  denominator  of  equation  (7)  is 


_ (0,0,0) _ 

k^Jmlm^yv]  +  mlmlv\  +  mimlu] 

\  X  y  s  X  ^  s  y  ^  s 

^^(0,0,0) 

4^  ^mlm^yCO^ {(p^)  +  mlml  sin"(^Jsin‘(0J  +  rrP'ynl  sin^(^Jcos"(^J 


(56) 


Then,  from  equation  (7),  the  DF  for  a  baffled  volumetric  aperture  is  approximated  as 

DF  =  sin'(0Jsin'(0J  +  sin^(^!i^)cos’(0^) ,  (57) 

At 

where  it  is  assumed  that  m^,  m^,  and  m,  are  all  significantly  greater  than  a  wavelength. 

The  product  above  may  be  considered  the  RMS  effective  area  of  the  collapsed 

weighting  in  the  xy-plane,  that  is, 

1/2 

(58) 


1 

7  r 

XT'  v/^(x)dx  J 

^y'v/yiy)dy 

j_w,(x)iir 

24 


where  =  J  J  w{x,y,z^dydz  diVid.  =  ^  J  w(x,>',z)<ix:rfe.  Notice  that  the  term 

cos(^^),  which  is  multiplied  by  the  product  m^m^.  in  equation  (57),  is  the  cosine  of  the  angle 

between  the  steering  direction  and  the  normal  to  the  xy-plane  (consider  a  unit  steering  vector 
e  =c^cos(0Jsin(0J  +  c^,sin(0Jsin(0J  +  e,cos(0J;then  e  e^  =cos(0j)).  Thus,  the  RMS  area 

of  the  collapsed  weighting  in  the  xy-plane,  as  seen  from  the  steering  direction,  is 

=m,mv|cos(0j|.  (59) 


In  a  similar  manner,  the  RMS  areas  of  the  collapsed  weightings  in  the  other  planes  are 

A,,  =m;^wjsin(0,)sin((/)j|,  (60) 

A,^  =m,.mjcos(0^)sin(0j|.  (61) 


Hence,  the  baffled  volumetric  directivity  approximation,  equation  (57),  returns  to  the 
familiar  form 


DF  = 


16;r~A^ 

A' 


16;r‘ 

~Y' 


(4  +  4  +  ^ 


2 


(62) 


where  A^  is  the  RMS  effective  aperture  area  for  the  given  steering  direction 

As  an  example,  consider  a  continuous  flat  weighting  function,  w(x,  y,  z)  =  constant,  for  a 
cubic  volumetric  aperture  |x|  |y|  <  LJl,  and  |z|  <  LJ2 .  Then,  from  equation  (53a), 


(63) 


The  effective  RMS  area  is  given  as 


4=^(L^Z^,cos^(^,)  +  L^L^sin^(0,)sin^((^,)  +  Z54cos^(0,)sin^(f)f',  (64) 
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and  the  approximate  directivity  of  the  baffled  cubic  aperture  with  flat  weighting  becomes 


DF  = 


cos^  ( (^  J  +  L^.  sin  ^  )  sin  ■  ( )  +  L^.  ll  cos"  ( 0^ )  sin " 


As  expected,  for  L,  =  0  and  <p,  =  0,  equation  (65)  reduces  immediately  to  twice  the  value  of 
equation  (39),  the  approximate  directivity  of  a  uniformly  weighted  nondirectional  planar  aperture  in 
the  free  field  (not  baffled). 


Equation  (57)  can  be  extended  to  include  discrete,  rather  than  continuous,  array  elements 
using  an  impulsive  weighting  function;  that  is. 


T 

w(x,y,z)  =  2^  w,  d(x  -  X, )  (5(>'  -  y, )  6(z  -  z, ) 
/=1 


(66) 


The  general  case  is  found  directly  by  substituting  equation  (66)  into  the  integrals  in  equation  (53a) 
and  evaluating  the  Gaussian  approximation  for  the  volumetric  DF  given  by  equation  (57). 

The  approximation  for  the  DF  of  a  volumetric  aperture  was  derived  by  using  a  Gaussian 
approximation  for  the  magnitude-squared  response.  The  result,  equation  (57),  uses  the  RMS 
lengths  of  the  collapsed  weightings  in  the  various  dimensions.  As  shown  in  appendix  A, 
alternative  approaches  for  linear  and  planar  apertures,  which  used  the  Parseval  relation,  yielded 
approximations  to  the  DF  that  depended  on  effective  lengths  or  areas,  respectively,  of  the 
weightings.  The  analogous  derivation  for  a  volumetric  aperture  or  array  has  not  been  discovered. 
Nevertheless,  by  extension,  the  following  conjectured  form  for  the  directivity  of  a  volumetric  array 
with  equal  element  spacing  d  in  all  three  dimensions,  based  upon  the  earlier  results  for  the  linear 
and  planar  arrays,  is  offered  (see  appendix  B). 

If  {w^.Xn,fn,k)}  are  the  three-dimensional  weights  applied  to  the  equispaced  volumetric 
array,  the  three  collapsed  (or  projected)  two-dimensional  weight  structures  are  defined  as 

k  m  n 
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The  three  effective  numbers  of  weights  are  then  defined  according  to 


■2 

2 

-  nk 

N=- 

_  nik 

nk 


mk 


(68) 


The  conjectured  directivity  at  operating  frequency/for  a  volumetric  array  is  then  given  by 


DF  =  4n^N,  =  +  Nlvl  + 

^  Jo 


\'A 


(69) 


Here,/^  =  c  /  {2d)  is  the  design  frequency,  where  the  common  element  spacing  J  is  a  half¬ 
wavelength.  The  steering  parameters  are  =  sin(0,)cos(ej,  =  sin(0^)sin(ej,  and 
=COS(0J. 
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3.  COMPARISON  WITH  EXACT  RESULTS 


In  this  section,  comparisons  are  made  between  the  exact  DF  (computed  by  numerical 
integration)  and  an  approximate  DF,  for  a  variety  of  weighted  linear,  planar,  and  volumetric  arrays. 
In  most  cases,  Hanning  weighting  is  employed,  and  the  array  elements  are  nondirectional. 

Section  2  gave  Gaussian  approximations  for  the  DFs  of  various  arrays.  However,  the 
Gaussian  approach  only  matches  the  mainlobe  region  of  the  array  power  response  and  takes  no 
account  of  the  sidelobes,  which  may  be  very  significant.  An  alternative  approach,  specifically  for 
arrays,  that  does  take  into  account  both  the  mainlobe  and  the  sidelobes  is  presented  in  appendix  B. 
This  latter  approximate  procedure  has  been  found  to  be  more  accurate  and  is  used  for  all  the 
numerical  examples  in  this  section.  As  will  be  seen,  the  approximations  are  quite  accurate,  even 
for  arrays  with  relatively  few  elements  and  steered  near  endfire.  With  increasing  array  size,  the 
approximations  become  more  accurate,  which  is  just  when  approximations  are  most  beneficial;  that 
is,  numerical  integration  can  become  time  consuming  and  burdensome  for  large  arrays. 

LINEAR  ARRAY 

Figure  6  shows  that  for  the  equispaced  line  array  with  spacing  d  and  steering  parameter 
=  cos  =  0.2,  the  approximation  for  the  DF  is  valid  beyond/=/^;  in  fact,  the  fit  is  good  up  to 
/=  1.6/^.  This  result  corresponds  to  element  spacing  J  =  1.6  XJl,  that  is,  greater  then  half¬ 
wavelength  spacing.  This  behavior  occurs  because  smaller  values  of  Im^I  keep  the  closest  grating 
lobe  of  the  array  response  at  greater  distances  from  the  visible  region,  thereby  allowing  the 
approximation  to  have  wider  applicability.  In  fact,  the  approximation  is  valid  for  <i<  2/(1  +  Im^I), 
which  corresponds  tof/f^  <  2/(1  +Im^I).  For  =  0.2,  the  latter  bound  isf/f^  <  1.67,  which  is  in 
good  agreement  with  figure  6.  For  the  equispaced  line  array  with  =  0.8,  the  bound  is 
f/f^  <1.1,  which  agrees  well  with  figure  7,  especially  for  the  larger  values  of  N. 
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Figure  6.  DF  for  an  Equispaced  Line  Array  with  Hanning  Weighting, 


Figure  7.  DF  for  an  Equispaced  Line  Array  with  Hanning  Weighting, 
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An  example  of  a  40-element  line  array  with  directional  amplitude  response  |  sin(0)sin(^)  | , 
that  is,  power  v  =  1 ,  is  given  in  figure  8.  Because  both  the  exact  result  and  the  approximation  for 
the  DF  involve  steering  angle  6^  only  through  the  multiplicative  term  sin‘(0J,  this  common 
quantity  is  factored  off;  this  leads  to  a  plot  of  DFy^sin‘(0^),  which  depends  only  on  the  remaining 

polar  steering  angle  0^. 

For  polar  steering  angle  0^,  near  0  or  ;r(endfire  steering),  the  parameter  =cos0^  is  near 
±1 .  The  approximate  DF  is  not  as  accurate  in  this  case  for  two  reasons.  The  function 
varies  fastest  at  the  edges  of  the  integration  range  (-1,  1)  on  u,  thereby  making  its  extraction  from 
the  integral  as  (1  -  ulY  less  accurate.  Also,  the  sidelobes  of  |pF(A:w- can  extend 

significantly  beyond  m  =  1  if  is  near  +1  (beyond  u  =  -\  if  is  near  -1);  these  sidelobes  would 
not  be  encountered  in  the  exact  integral  for  V  over  (-1, 1).  However,  the  aproximating  integral  for 
V  extends  over  all  the  sidelobes  on  both  sides  of  the  mainlobe,  halfway  up  to  the  next  grating  lobe. 

This  lack  of  accuracy  for  0^  near  0  and  K  is  illustrated  in  figures  8  and  9,  where  the 
discrepancy  of  the  exact  DF  from  the  approximate  DF  is  significant  even  when  ///„  <  1 ,  that  is, 
d  <  A/2 .  The  errors  are  less  for  Hanning  weighting  than  for  flat  weighting  because  Hanning 
sidelobes  decay  to  low  levels  faster,  tending  to  minimize  the  two  effects  discussed  above.  Of 
course,  for  both  weightings,  the  approximate  DF  only  applies  for  ///„  <  2/(l-t- 1  |)  or 

dj X<  l/(l-t- 1  Wj  I)  =  l/(l+ 1  cos(0^)  l),  at  which  point  grating  lobes  begin  to  manifest  themselves. 
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►XIMATION 


'h  Hanning  Weigh 


APPROXIMATION 


■71 


The  next  example  consists  of  a  curved  line  array,  specifically  an  equispaced  semicircular 
array  centered  at  the  origin  in  the  xy-plane  with  its  middle  element  located  on  the  jc-axis;  that  is,  the 
angular  locations  of  the  N  elements  are 


for  1  <  «  <  iV .  The  elements  are  all  located  at  radius  r  and  employ  flat  weighting.  Although  this 
case  might  be  considered  a  (sparse)  planar  array,  its  DF  behaves  like  that  of  a  line  array,  being 
linearly  proportional  to  the  effective  curved  line  length,  namely. 


2  N  ^  N  r 

—Ttr - =  iTt - . 

A  iV-1  N-\  A 


(The  factor  N/{N  -1)  represents  the  effective  extension  of  the  edge  elements,  when  the  weighting 
is  flat.)  The  approximate  DF  is  plotted  as  a  dashed  line  in  figure  10  and  agrees  very  well  with  the 
exact  results  for  several  polar  steering  angles.  The  tic  mark  labeled  “H”  indicates  where  the 
distance  between  adjacent  elements  is  a  half-wavelength;  the  approximate  DF  is  excellent  for 
smaller  spacings,  despite  the  fact  that  the  sidelobes  of  this  particular  unweighted  array  are  rather 
poor  (~  -8  dB).  The  breakpoint  between  the  approximation  and  exact  results  is  located 
approximately  where  the  spacing/wavelength  ratio  is  2/(\+ 1  |)  =  2/(1  -i-sin0^)  for  this  example 

where  azimuth  steering  angle  =  0. 

The  only  difference  in  figure  1 1  is  that  azimuthal  steering  angle  6^  is  changed  to  k/I,  that 

is,  along  the  y-axis.  Again,  the  approximate  DF  above  is  excellent  when  the  spacing/wavelength 
ratio  is  less  than  2/(l  +  max{|  |,|  |})  =  2/(1  + sin 0,)  for  e,  =  nl2. 
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Figure  11.  DF  for  an  Equispaced  Semicircle  Atray  with  Flat  Weighting, 

N  =  48,  d  =  K  /2 


PLANAR  ARRAY 


For  a  baffled  planar  array  with  elements  spaced  at  distance  d  in  both  the  x  and  y 
dimensions,  the  upper  bound  on  operating  frequency/,  for  which  the  DF  approximation  is  valid,  is 

/  2 

<  - 

/o  l^i^aY^+b 

where  a  =  min(|  |,|  |),  b  =  max(|  w,  |,|  |) .  This  result  agrees  well  with  the  baffled  planar  array 

results  in  figures  12  and  13.  A  more  accurate  bound  is  given  in  appendix  B. 


Figure  12.  DF  for  an  Equispaced  Baffled  Planar  Array  with  Multiplicative 
Hanning  Weighting,  d^  =  d^,  u^  =  0.5,  =  0.6 
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Figure  13.  DF  for  an  Eguispaced  Baffled  Planar  Array  with  Multiplicative 
Hanning  Weighting,  d^  =  d^,  u^  =  0.2,  =  0.3 

VOLUMETRIC  ARRAY 

An  example  of  the  application  of  the  DF  approximation  (equation  (69))  is  given  in  figure 
14  for  an  equispaced  volumetric  array  with  multiplicative  Hanning  weighting  in  the  three 
dimensions.  The  numbers  of  elements  are  =  10,  tV^.  =  12,  AC  =  15,  and  the  steering  direction  is 
=  0,  Vs  =  0,  =  1 ;  that  is,  toward  the  North  Pole. 

Results  for  an  unbaffled,  as  well  as  a  baffled,  volumetric  array  are  presented  in  figure  14 
(solid  curves),  along  with  the  directivity  approximation  (dashed  curve).  The  agreement  between 
the  exact  baffled  DF  and  the  conjectured  DF  is  excellent  up  to  approximately  f  / fo  -  1-8,  whereas 
the  exact  unbaffled  DF  differs  significantly  in  the  neighborhood  of//yC  =  1  and  ioxf/fo>  1.7. 
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Figure  14.  DF  for  an  Equispaced  Volumetric  Array  with  Multiplicative 
Hanning  Weighting,  u^  =  0,  Vs  =  0,  Ws  =  1 


These  differences  occur  because  of  the  locations  of  the  aliasing  lobes  of  the  unbaffled 
volumetric  array.  For  the  present  equispaced  array  with  common  separation  distance  d  in  three 
dimensions  and  separable  weighting,  the  array  function  is  given  by  the  product  of  three  functions 
of  the  form 


'^w^„tx^inkd{u-u)j  , 

n=I 


kd  =  n-- 
fo 


(70) 


They  and  z  functions  are  similar,  with  replacements  v- and  w  -  respectively,  inside  the 
exponential.  This  function  of  u  has  period  2k  !  {kd)  =  Ift  f  with  the  desired  mainlobe  at  w  =  u^. 
The  remaining  aliasing  lobes  at  ±  m2fl /can  become  major  grating  lobe  contributors  under 

conditions  to  be  developed  below. 

If  the  double  integral  on  6, 0  in  equation  (7)  is  interpreted  as  a  triple  integral  in  u,  v,  w 
space,  the  integration  takes  place  only  on  the  unit  sphere.  Therefore,  significant  contributions  to 
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he  integral  can  take  place  only  when  the  three  array  functions  have  combinations  of  mainlobes 
and  at  least  one  aliasing  lobe  that  intersect  on  the  sphere.  For  example,  one  such  case  is 


->  7 

U'  +V'  + 


f 

w  -  2^ 
^  / 


=  1; 


that  is, 


/ 


(71) 


More  generally,  the  frequencies  at  which  a  significant  contribution  occurs  to  the  double  integral  in 
equation  (7),  and  only  one  aliasing  lobe  is  involved,  are 

/I  1  1  ni\ 

—  =  i — I  or  —  or -  (/2) 

fo  ^5  '^V 


When  two  aliasing  lobes  are  involved,  the  interfering  frequencies  are 

/  2  2  2 

—  _ _  _  r\-r  -  . 


or  - 

or  - 

+  Vs 

W,  + 

(73) 


Finally,  when  three  aliasing  lobes  are  involved  (no  mainlobes),  the  interfering  frequencies  are 

/  3 

fo  +  Vs  + 


In  general,  the  lowest  frequency  at  which  grating  lobes  contribute,  for  an  unbaffled  equispaced 
volumetric  array,  is 

Ism.  = _ I _  (75) 

fo  max{|wj,|i/j,|wj} 

For  example,  in  figure  14,  with  u^,  v^,  =  0,  0, 1,  the  problematical  frequencies  for  an 

unbaffled  array  are f  / /„  =  I  /  I  and 2  /  (w5+h'j)  =  2,  at  least  for  0  <//^  <  2.  The  lower 
frequency  corresponds  to  a  lower-hemisphere  response  at  w  =  w^-2f^l f  =  l-  2  =  -l.  Because 
the  baffled  array  cannot  respond  in  the  lower  hemisphere,  it  has  no  major  response  at//^  =  1, 
thereby  maintaining  the  larger  value  of  directivity,  in  agreement  with  the  approximate  DF  (which 
ignores  grating  lobes). 
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The  width  of  the  discrepancy  between  exact  and  approximate  DFs  is  related  to  the 
beamwidth  of  the  array  in  the  angular  dimensions.  Larger  numbers  of  elements  will  lead  to 
narrower  beamwidths  and  therefore,  narrower  regions  of  discrepancy  between  the  unbaffled  array 
DF  and  the  approximate  DF. 

The  top  example  in  figure  15  uses  an  unbaffled  volumetric  array  with  N„  Ny,  =  20,  24, 
30  elements,  respectively,  and  steering  direction  v^,  =  0.6,  0.4,  0.693.  The  lowest 

interfering  frequency  occurs  at//^  =  1  /  =  1 .44,  due  to  the  aliasing  lobe  atw  =  Ws-2fo/f= 

-0.693.  This  is  a  lower-hemisphere  response. 

The  next  lowest  frequency  is  // /„  =  2/(w^  -i-  wj  =  1.55.  The  contributing  lobes  are  at 
u  =  u^-2f^l  f  -  -0.693,  v=  Vs  =  0.4,  and  'w-w^-2f^/ f  =  -0.6.  Notice  that  the  sum  of 
squares  of  these  three  locations  is  1,  meaning  on  the  unit  sphere.  This  is  also  a  lower-hemisphere 
response. 

The  next  frequency  is  ///„  =  l/  =  1.67.  The  contributing  lobes  are  at 
u=-u^-2f^l  f  =  -0.6,  v=  Vs  =  0.4,  and  w  =  =  0.693.  This  is  an  upper-hemisphere  response. 

(For  a  baffled  array,  this  would  be  the  lowest  frequency  leading  to  an  aliased  response.) 

Finally,  there  are  two  other  responses  that  blend  together,  namely  flfo"=2l{Us+Vs+ 

=  1.77  and  2  /  (v^  +  =  1 .83.  The  aliasing  lobe  locations  are  at  w  =  -  2^  //=  -0.53, 

V  =  Vs-2fol /=  -0.73,  w  =  Ws-2fol  f=  -0.437  for  the  first  case,  and  at  w  =  =  0.6,  v  =  v^  -  2/^  // 

=  -0.693,  w  =  Ws-2fol /=  -0.4  for  the  second  case.  Both  of  these  cases  are  lower-hemisphere 
responses. 

The  bottom  imbaffled  example  in  figure  15  uses  eight  times  fewer  elements,  namely  1800 
elements  instead  of  14,400.  The  decreased  resolution  capability  causes  several  of  the  aliased 
responses  to  overlap,  making  it  impossible  to  identify  individual  aliased  contributions. 
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Figure  15.  DF  for  an  Equispaced  Unbaffled  Volumetric  Array  with  Multiplicative 
Hanning  Weighting,  u^  =  0.6,  =  0.4,  =  0.693 

APPLICABILITY  TO  VOLUMETRIC  ARRAYS,  BAFFLED  AND  UNBAFFLED 

Although  there  is  only  one  form  for  the  approximate  DF  of  a  volumetric  array,  this  form 
has  different  ranges  of  applicability,  depending  on  whether  baffling  is  used.  The  approximate  DF 
takes  no  account  of  grating  lobes,  whereas  actual  arrays  are  always  subject  to  the  effects  of  grating 
lobes.  This  is  evident  in  figure  14,  where  the  agreement  of  the  approximate  DF  with  the  exact 
baffled  DF  is  excellent  up  to  ///„  =  1 .8,  whereas  the  agreement  with  the  exact  unbaffled  DF 
holds  only  up  to  ///,  =  0.9 .  Baffling  suppresses  all  the  lower-hemisphere  grating  lobes;  the  dip 
in  the  unbaffled  DF  near  ///„  =  1  is  due  to  a  lower-hemisphere  grating  lobe  which  is  responding 
to  undesired  noise  arrivals.  However,  it  should  be  noted  that  for  frequencies  1.1<  f  j  f^<\.l , 
the  approximate  DF  again  agrees  with  the  exact  unbaffled  DF;  this  occurs  because  the  particular 
grating  lobe  in  figure  14  only  contributes  in  the  neighborhood  of  ///„  =  1. 
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The  general  frequencies  at  which  an  equispaced  unbaffled  volumetric  array  is  subject  to  a 
grating  lobe  are 

fc  +  + 

where  I,  J,  and  K  are  positive  and  negative  integers,  not  all  equal  to  zero.  The  lowest  frequency 
was  given  earlier  in  equation  (75). 

For  an  equispaced  baffled  volumetric  array,  not  all  the  frequencies  in  equation  (76)  lead  to 
grating  lobes.  To  determine  which  frequencies  are  relevant,  a  search  over  all  the  small  integer 
values  for  7, 7,  and  K,  such  as  0,  ±1,  ±2,  must  first  be  conducted.  Then,  those  cases 
corresponding  to  negative  values  of  +  Klfjf  must  be  discarded,  because  they  correspond  to 

lower-hemisphere  grating  lobes,  which  will  not  be  active  in  a  baffled  application.  The  remaining 
cases  constitute  valid  grating  lobes;  the  lowest  remaining  frequency  will  be  the  upper  bound  on 
applicability  of  the  approximate  DF  to  an  equispaced  baffled  volumetric  array. 
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4.  SUMMARY 


Estimating  the  directivity  of  an  array  remains  a  fundamental  and  important  aspect  of  array 
design.  These  estimates  are  used  in  many  applications:  to  determine  frequency  regions  in  which 
an  array  will  be  ambient  or  self-noise  limited,  to  determine  detection  ranges,  and  as  a  comparison 
baseline  for  various  changes  in  array  parameters  (such  as  geometry,  element  shading  coefficients, 
and  noise  distribution).  It  is  necessary,  particularly  when  one  is  trying  to  examine  the  merits  of  a 
given  set  of  conditions  over  another,  to  have  accurate  estimates  for  directivity.  The  approximations 
presented  in  this  report  for  the  DFs  of  linear,  planar,  and  volumetric  apertures  and  arrays  provide 
quick  and  reasonably  accurate  directivity  estimates  without  having  to  resort  to  numerical 
integration  and  time-consuming  computer  calculations. 

In  figure  2,  it  was  shown  that  a  commonly  used  approximation  for  the  directivity  of  a  line 
array  is  valid  only  with  omnidirectional  elements  and  at  the  design  frequency  /„,  that  is,  at  a  single 
frequency.  At  other  frequencies,  the  DF  becomes  a  function  of  steering  angle,  and  the  standard 
approximation  is  inaccurate.  The  DF  approximations  derived  here  are  more  general  and  allow  for 
variable  steering  angles  and  frequencies  of  operation. 

Another  useful  result  of  this  study  is  illustrated  in  figures  3  and  5,  which  show  that  highly 
directional  array  elements,  such  as  elements  having  an  angular  response  proportional  to  |  cos(0)  |'', 
can  lead  to  significant  decreases  in  array  directivity  at  moderate  steering  angles  from  broadside. 
Therefore,  it  may  be  surmised  that  element  directivity,  for  most  steering  angles,  reduces  directivity, 
although  one  must  bear  in  mind  that  these  calculations  assume  a  spherically  isotropic  noise  field. 
Directional  elements,  in  general,  improve  array  gain  for  steering  angles  near  broadside. 

The  directivity  factor  of  a  linear  aperture  or  array  depends  on  the  first  power  of  the  effective 
length  of  the  line.  This  results  from  the  inability  of  a  (vertical)  line  to  resolve  in  the  azimuthal 
direction;  that  is,  it  has  a  conical  response  region,  centered  on  the  line,  within  which  all  arrivals  are 
treated  alike.  This  inherent  ambiguity  of  a  line  keeps  the  DF  at  moderate  values.  However,  this 
lack  of  resolution  can  be  partially  remedied  by  using  directional  elements  in  the  line  aperture  or 
array,  leading  to  a  modest  increase  in  the  DF. 

For  a  planar  or  volumetric  aperture  or  array,  the  situation  is  markedly  improved,  because 
simultaneous  resolution  in  both  the  azimuthal  and  polar  angles  is  possible.  In  fact,  when  the 
effective  area  is  large  compared  to  the  squared  wavelength,  rather  good  resolution  can  be  achieved; 
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this  fine  resolution  overrides  any  modest  directivity  that  individual  elements  may  have,  making  the 
overall  DF  virtually  independent  of  the  individual  element  properties. 

Baffling  a  planar  aperture  or  array  doubles  the  DF  by  virtue  of  eliminating  any  arrivals  in 
the  undesired  hemisphere.  On  the  other  hand,  a  volumetric  aperture  or  array  has  the  inherent 
capability  of  being  able  to  reject  arrivals  in  an  undesired  hemisphere,  and  therefore  need  not 
employ  baffling.  The  DF  in  the  volumetric  case  is  given  by  the  same  formula  as  for  a  baffled 
planar  array,  but  where  the  effective  area  must  be  interpreted  as  that  seen  from  the  steering 
direction.  In  addition  to  steering  direction,  the  effective  area  will  incorporate  the  three-dimensional 
physical  extent  and  detailed  weighting  of  the  volumetric  aperture  or  array. 

The  DF  approximations  for  linear,  planar,  and  volumetric  apertures  and  arrays  are 
summarized  in  appendix  C.  Results  obtained  by  the  use  of  the  Gaussian  approximation,  as  well  as 
ParsevaFs  relation  and  a  modification  for  arrays,  are  listed. 
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APPENDIX  A 

ALTERNATIVE  INTEGRAL  APPROXIMATIONS 

LINEAR  APERTURE 

From  equation  (20),  the  amplitude  response  of  a  continuously  weighted  directional  line 
aperture  is  given  as 

OO 

A(d,(t>)  =  |sin(0)sin((^)|'^  |  w(z)exp|zfe(cos(0)  -  cos(  .  (A-1) 


The  value  v  =  0  corresponds  to  an  omnidirectional  element  angular  response.  The  directivity 
factor  is  again  expressed  as 


DF  = 


4n\A{d„<t>,f 


iTtK 

Jj|A(0,0)f  sm{^)d<pde 


0  0 


(A-2) 


Then  the  integral  which  must  be  approximated  is,  as  in  equation  (21), 


Ik  K 

=  J |sin^(0)|  d9^ |sin^(0)|  sin(0) | w(z)exp|/fe(cos(^) -  cos(<p^y)^di 


(A-3) 


Expressing  the  magnitude-squared  as  a  double  integral  on  z  and  ^  and  letting  u  =  cos(0)and 
s  =  z-^,  where  ^  is  a  dummy  variable  of  integration,  equation  (A-3)  becomes 


I  =  h(v)  J  J w(z)w(z  -  5')expj-zfoM^|J (1  -  u^Y  exp^iksu^dudzds , 


(A-4) 


where  h{v)  is  defined  by  equation  (11).  The  inner  integral  over  u  in  equation  (A-4)  may  be 
evaluated  exactlyi^-  P-321]  yielding 


OO  OO  _ 

I  =  h(v)  J  J H'(z)w(z  -  s)cxp^-iksu^^^r(v  -t- 1)  ‘  J^^p_{ks)dzds.  (A-5) 
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By  using  the  relation  h{v)  -  r(v  +  l/2)/r(v  + 1)  in  equation  (1 1)  and  the  weighting 
autocorrelation  function  defined  as 


=  I  H’(z)w(z  -  s)dz,  (A-6) 


equation  (A-5)  can  be  expressed  as 


7  =  2;rr(v  + 1/2)  j<I>^(5)exp|-/focos(^j| 


ks 


(A-7) 


However,  because  5)  varies  slowly  with  5  relative  to  the  remaining  functions  in  equation 
(A-7),  the  following  approximation  obtains: 


7  =  2""'/';^r(t;  +  l/2)<D,,(0)J 


cos(focos(0J)J,^,,2(^^) 


{ks) 


v+1/2 


■ds. 


(A-8) 


using  the  even  character  of  entire  function  J^{z)fz’^  in  z.  Finally,  use  of  Erdelyf  yields  the  result 


/  =  4^3/2  r(v  + 1/2)  ^in2v(^  )  ^ 

ytr(v  +  l) 


The  approximate  DF  follows  from  equation  (A-2)  according  to 


(A-9) 


24  V^nY_^g-^2v^g  ^ 
X  r(v+i/2) 


where  L  is  the  effective  length  of  the  weighting. 


(A- 10) 


L  = 


j  w{z)dz 


e  ^ 


(A-ll) 


jw^(z)dz 
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In  the  special  case  of  a  nondirectional  aperture,  v  =  0 ,  equation  (A- 10)  simplifies  to  the  standard 
result 


DF  = 


(A-12) 


which  is  independent  of  the  steering  direction. 


More  generally,  for  arbitrary  element  directional  response  g(0, 0)  in  equation  (A-1)  instead 
of  I  sin(6)sin(0)  f,  the  approximate  DF  in  equation  (A-10)  is  replaced  by 


(A-13) 


A  more  direct  derivation  of  the  result  in  equation  (A-9)  may  be  obtained  by  use  of 
Parseval’s  relation:  from  equations  (22)  and  (23),  the  integral  on  u  is,  for  |  |=|  cos(0^)  |<  1, 

j|PF(A:(w-Mj)|  du^\y[k{u-uS\  du  =  l\w\z)dz  =  (A-14) 

-1  -oo  -oo 

Then,  equation  (22)  yields 

/  =  ;t(v)sin^''(^JA<DJO),  (A-15) 

which  is  identical  to  equation  (A-9).  The  same  approximation  for  the  DF,  as  given  by  equations 
(A-10)  and  (A-1 1),  results. 
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PLANAR  APERTURE 


An  alternative  to  expression  (38)  for  the  DF  of  an  unbaffled  planar  aperture  is  possible 
through  the  use  of  Parseval’s  relation  in  two  dimensions.  Namely,  from  equations  (32)  and  (35), 
there  follows  immediately 


OO 

7  s  2(cos(0j)’’'^  'A"Jj  Mr{x,y)dxdy 


(A- 16) 


for  the  denominator  of  equation  (7).  The  approximate  DF  then  readily  follows  as 
DF  =  2;r^cos(^,), 


(A-17) 


where  the  effective  area  of  weighting  w(x,y)  is  defined  as 


4  = 


n2 


Jj  w{x,y)dxdy 


Jj  w^{x,y)dxdy 


(A-18) 


(For  uniform  weighting  over  physical  lengths  and  L,.,  this  definition  reduces  to  =  L^L^,  as 
noted  under  equation  (39)). 
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APPENDIX  B 

ALTERNATIVE  APPROXIMATE  DIRECTIVITY  FACTOR 
FOR  EQUISPACED  ARRAYS 


In  the  main  text,  an  approximate  DF  was  derived  for  linear  and  planar  apertures  by 
matching  the  mainlobe  region  of  the  beam  pattern  power  response  with  a  Gaussian  function.  This 
approximate  result  was  then  extended  to  discrete  linear  and  planar  arrays  by  replacing  the 
continuous  weighting  functions  by  sets  of  weighted  impulses.  The  end  results  employed  RMS 
measures  of  the  extent  of  the  aperture  or  array. 

The  major  limitation  of  the  Gaussian  approximation  is  that  it  takes  no  account  of  the  actual 
sidelobe  levels  of  the  beam  pattern;  that  is,  the  Gaussian  approximation  matches  only  the  three 
lowest-order  derivatives  of  the  response  in  the  mainlobe  region.  This  procedure  could  be  extended 
to  higher  derivatives  by  using  a  more  general  fit,  such  as  exp(-ax:^)(^  +Bx^y,  however,  this 
modification  would  still  not  take  account  of  the  distant  sidelobes. 

This  appendix  presents  an  alternative  approximation  for  the  DF  of  equispaced  linear  and 
planar  arrays  that  does  account  for  the  sidelobes.  (It  does  not,  however,  account  for  the  grating 
lobes.)  Then,  the  result  is  extended  to  equispaced  volumetric  arrays  by  means  of  a  conjectured 
form.  This  new  procedure  can  be  thought  of  as  an  extension  of  Parseval’s  relation  (which  is 
applicable  to  apertures  with  continuous  weighting)  to  discrete  arrays  with  impulsive  weighting. 

The  end  results  will  employ  effective  measures  of  the  array  extent. 

EQUIPSACED  LINEAR  ARRAY 

In  equation  (22),  the  double  integral  over  angles  ^,6  ,for  the  noise  output  power  of  a  linear 
aperture,  was  simplified  to 

1  2 

V  =  2K^^{ku-ku^du ,  W(v)  =  jw(z)exp(ivz)dz ,  (B-1) 

-1 

for  continuous  weighting  w(z)  and  v  =  0,  where  steering  parameter  =  cos(p^.  For  an 
equispaced  weighted  line  array  on  the  z-axis,  the  weighting  takes  the  special  form 

w{z)  =  ^w„6iz-nd),  (B-2) 


B-1 


leading  to  window  function 


W(v)  =  '^w^  expiindv) 

n 


(B-3) 


for  all  V,  which  has  period  Injd  in  v. 


The  integral  in  equation  (B-1)  now  takes  the  form 
1 

V^2k\ 


w„exp| 


f .  iTtd 


in- 


V 


A 


du . 


(B-4) 


The  integrand  has  a  desired  peak  at  u  =  u^,  along  with  undesired  periodic  grating  lobes  separated 
from  by  multiples  of  Xfd.  If  element  spacing  d  is  less  than  a  half-wavelength,  then  the  period 

P=Xld  >2,  and  all  the  grating  lobes  of  ]^{ku-ku^'f  lie  outside  of  the  (-1,1)  interval  in  u. 
Therefore,  if  the  array  length  is  much  greater  than  a  wavelength,  and  if  the  sidelobes  are  low  level, 
the  integration  interval  in  equation  (B-4)  can  be  modified,  leading  to  approximation 


Pf2 


-Pfl 


(  .  2mi  \ 
exp  m—u 


du  =  2;r— V 

d^  " 


(B-5) 


(This  last  quantity  is  the  discrete  analogue  to  Parseval’s  relation  for  apertures,  which  yields  the 
integral  of  the  squared  weighting.)  Use  of  relations  (B-3)  and  (B-5)  now  yields  the  DF 
approximation  for  a  linear  equispaced  array  in  the  form 


DF  = 


AkW^{0)  _  2d 
V  ~  X  ' 


( 


n 


(B-6) 


The  quantity  is  the  effective  number  of  weights  in  real  set  j:  (For  flat  weights,  is 
simply  the  number  of  elements.)  If  quantity  dK  is  interpreted  as  the  effective  length  of  the 
discrete  array,  relation  (B-6)  for  the  approximate  DF  takes  the  same  form  as  for  a  linear  aperture, 
namely,  2  LgfX. 
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For  elements  with  amplitude  directionality  |sin0  sin ^|,  relation  (B-6)  is  somewhat 
modified,  to  the  extent  that  the  factor  of  2  in  the  DF  is  replaced  by  4  sin  6^ .  The  factor  is 
unchanged. 

Although  the  derivation  above  presumed  dIX  <  1  /  2 ,  result  (B-6)  is  actually  valid  for  the 
broader  region  djX  <  l/(l-i- 1  I).  This  is  due  to  the  fact  that  the  nearest  grating  lobe  does  not  yet 

affect  the  noise  integral  for  this  region  of  values  of  djX. 

EQUISPACED  PLANAR  ARRAY 

In  equation  (34),  the  double  integral  over  angles  ^,0 ,  for  the  noise  output  power  of  a 
baffled  planar  aperture  in  the  xy-plane,  was  simplified  (for  v  =  0)  to 

1 

V  =  — Up  ^  ’  (B-7) 

c, 

where  C]  is  a  unit  circle  centered  at  the  origin,  steering  parameter  =  cos<p^  >  0,  and  window 

W{u,  v)  =  jj  w(x,>^)exp(/ux  -I-  ivy)dxdy  (B-8) 

for  all  u,v,  which  depends  on  continuous  two-dimensional  weighting  w(x,  y).  For  an  equispaced 
weighted  planar  array,  the  weighting  takes  the  special  form 

w{x,y)  = 

nm 

leading  to  window  function 

W{u,v)  =  exp(mJ^M  +  irnd^v)  (B-10) 

nm 

for  all  u,v,  which  has  period  2k Jd^  in  u  and  period  Injdy  in  v. 

Substitution  of  equation  (B-10)  in  equation  (B-7)  yields 
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v  =  -!-jJ 


Ijld 

in^{u-u^)  +  im-:^{v-v^) 
A  A 


dudv. 


The  integrand  has  a  desired  peak  at  u  =  u^,v=v^,  along  with  undesired  periodic  grating  lobes 
separated  from  by  multiples  of  Kjd^  in  u,v,  respectively.  If  element  spacings 

d^  and  dy  are  both  less  than  a  half-wavelength,  then  the  periods  Xfdy.  and  AJdy  are  both  greater 

than  2,  and  all  the  grating  lobes  of  ^{ku-ku^  ,  kv-  kv^  jj^  lie  outside  of  the  unit  circle  Q  in  the 

MV-plane.  Therefore,  if  the  array  lengths  in  the  x-  and  y-directions  are  both  much  greater  than  a 
wavelength,  and  if  the  sidelobes  are  low  level,  the  integration  region  in  equation  (B-ll)  can  be 
modified,  leading  to  approximation 
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where  i?  is  a  rectangle  of  size  ±A  l{2d^)  by  ±  A  /(2<i^) ,  centered  at  the  origin  of  the  M\/-plane. 

(The  last  quantity  in  equation  (B-12)  is  the  discrete  analogue  to  Parseval’s  two-dimensional  relation 
for  apertures,  which  yields  the  integral  of  the  squared  weighting  >v(x,  y).)  Use  of  relations  (B-8) 
and  (B-12)  now  yields  the  DF  approximation  for  an  equispaced  baffled  planar  array  as 


Z)F  = 


4;rW"(0,0) 


=  COS0, 


1 4nd^d^. 
'  A' 


A7  =2^ _ 

X\-  2 

> 


(B-13) 


The  quantity  N^y  is  the  effective  number  of  weights  in  two-dimensional  real  set  (For  flat 

weights,  is  simply  the  total  number  of  elements  employed  in  the  xy-plane.)  If  the  quantity 
d,dyN^  is  interpreted  as  the  effective  area  of  the  discrete  planar  array,  relation  (B-13)  for  the 
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approximate  DF  of  a  baffled  planar  array  takes  the  same  form  as  for  a  baffled  planar  aperture, 
namely  |cos^^|4;r  /A^ .  An  alternative  expression  for  the  DF  of  a  baffled  planar  aperture  is 


DF  s4;r 


(B-14) 


where  Ag  =  |cos0^  \a^  is  the  effective  area  of  the  planar  aperture,  as  seen  from  the  steering 
direction  0^,(1)^. 


For  elements  with  moderate  amplitude  directionality,  relation  (B-13)  for  the  approximate 
DF  is  unchanged  because  a  two-dimensional  array  can  have  simultaneous  fine  resolution  ability  in 
both  the  azimuthal  and  polar  angles.  This  fine  resolution  dominates  the  relative  slowly  varying 
selectivity  of  the  individual  element  directionality,  effectively  factoring  it  out  of  the  double  integral 
on  0,0,  with  the  element  directionality  value  taken  at  the  steering  angle  6^,  05. 


Although  the  derivation  above  presumed  d^<  A  /2  and  dy  <  A/  2 ,  result  (B-13)  is 
actually  valid  for  the  broader  region  where 
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(B-15) 


for  all  integers  7  and  7  equal  to  -1,  0,  -t-1,  except  that  7,  7  =  0,  0  is  disallowed.  Restriction  (B-15) 
guarantees  that  none  of  the  grating  lobes  in  equation  (B-1 1)  lie  inside  the  unit  circle  Q;  this  is 

consistent  with  the  rectangular  replacement  used  in  approximation  (B-1 2),  which  guaranteed  the 
exclusion  of  all  grating  lobes  from  the  rectangular  integration  region. 
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EQUISPACED  VOLUMETRIC  ARRAY 


The  authors  have  been  unable  to  discover  how  to  extend  the  arguments  above  to  a 
volumetric  aperture  or  equispaced  array;  accordingly,  conjectured  forms  for  the  DF  of  volumetric 
apertures  and  arrays  are  constructed  based  on  the  earlier  forms  in  this  appendix  for  the  approximate 
DFs  in  the  planar  case. 

Whereas  a  planar  aperture  or  array  occupies  only  the  xy-plane,  a  volumetric  aperture  or 
array  occupies  the  third  dimension,  namely  z.  This  gives  the  steered  volumetric  amplitude 
response  the  capability  to  reject  responses  in  an  undesired  hemisphere,  whether  it  is  baffled  or  not. 
Therefore,  the  baffled  planar  DF  result  in  equation  (B-14)  should  furnish  the  model  form  for  the 
three-dimensional  cases  (whether  baffled  or  not),  namely 

(B-16) 


where  Ag  is  now  the  effective  area  of  the  volumetric  aperture  or  array,  as  seen  from  the  steering 
direction  6^(p^. 


In  particular,  guided  by  the  Gaussian  form  in  equations  (59)  through  (62),  there  follows, 
for  a  volumetric  aperture. 


=(' 


■ 


(B-17) 


where  the  projected  effective  areas  are  defined  in  terms  of  the  three-dimensional  aperture  weighting 
w(x,  y,  z)  according  to 


x\ 


\\'^^{x^y)dxdy 
jj  w%{x,y)dxdy 


(B-18) 


For  an  equispaced  volumetric  array  with  discrete  three-dimensional  weighting  w(n,m,k) 
and  spacings  d^,d^.,d,,  aproximate  DF  (B-16)  is  still  appropriate,  but  where  the  effective  area  is 

now  represented  as 
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(B-19) 


The  projected  effective  numbers  of  elements  here  are  given  in  terms  of  the  three-dimensional 
weighting  {yv{n,m,  *)}  as 

r  Y 

'  2  , - w^in,m)  =  yw{n,m,k),  etc.  (B-20) 


Although  there  is  only  one  form,  namely  equation  (B-16)  with  definitions  (B-19)  and 
(B-20),  for  the  approximate  DF  of  an  equispaced  volumetric  array,  this  form  has  different  ranges 
of  applicability,  depending  on  whether  baffling  is  used.  The  approximate  DF  takes  no  account  of 
grating  lobes,  whereas  actual  arrays  are  always  subject  to  the  effects  of  grating  lobes.  This  is 
evident  in  figure  14  where  the  agreement  of  the  approximate  DF  with  the  exact  baffled  DF  is 
excellent  up  to  f  jfo  =  1-8,  whereas  the  agreement  with  the  exact  unbaffled  DF  holds  only  up  to 
f  If  o  =  0-9-  This  is  due  to  the  fact  that  baffling  suppresses  all  the  lower-hemisphere  grating  lobes; 
the  dip  in  the  unbaffled  DF  near  /  ff^  =  1  is  due  to  a  lower-hemisphere  grating  lobe  that  is 

responding  to  undesired  noise  arrivals.  However,  it  should  be  noted  that  for  frequencies 
1 . 1  <  ///o  <  1 .7,  the  approximate  DF  again  agrees  with  the  exact  unbaffled  DF;  this  occurs 
because  this  particular  grating  lobe  in  figure  14  only  contributes  in  the  neighborhood  of  f  =  1. 

The  general  frequencies  at  which  an  equispaced  unbaffled  volumetric  array  is  subject  to  a 
grating  lobe  are 

/  1^+J^  +  K^ 

-r=i - h  (B-21) 

Jo  +v^J  +  w^K\ 

where  I,  J,  K  are  positive  and  negative  integers,  not  all  equal  to  zero.  The  lowest  problematical 
frequency  was  given  earlier  in  equation  (75). 

For  an  equispaced  baffled  volumetric  array,  not  all  the  frequencies  in  equation  (B-21)  lead 
to  grating  lobes.  To  determine  which  frequencies  are  relevant,  a  search  over  all  the  small  integer 
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values  for  I,  J,  K,  such  as  0,  ±1,  +2,  must  first  be  conducted.  Then,  those  cases  corresponding  to 
negative  values  of  +  K  2/^//  must  be  discarded,  because  they  correspond  to  lower- 

hemisphere  grating  lobes,  which  will  not  be  active  in  a  baffled  application.  The  remaining  cases 
constitute  valid  grating  lobes;  the  lowest  remaining  frequency  furnishes  the  upper  bound  on 
applicability  of  the  approximate  DF  to  an  equispaced  baffled  volumetric  array. 
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APPENDIX  C 

SUMMARY  OF  FORMULAS  FOR  APPROXIMATE  DIRECTIVITY  FACTOR 

The  approximate  relations  for  the  DF  are  listed  below.  In  each  case,  the  Gaussian  results 
are  presented  first;  then,  the  results  of  the  alternative  method  employing  Parseval’s  relation,  or  a 
modification  for  equispaced  arrays,  are  summarized. 

LINEAR  APERTURE 


m  ■ 


J  z^w(z)dz 
j*  w{z)dz 

J—oo 


Nl/! 


....  24  2ng‘(eM 

I  r2;r  ^ 

^  i  g\e,,i,,)de 


L  = 


J  w(z)dz 


Jw^(z)Jz 


LINEAR  ARRAY 


2d 

DF  =  —  N.. 
A  - 


(C-2) 


(C-3) 


(C-4) 


(C-5) 


(C-6) 


C-1 
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PLANAR  APERTURE  (UNBAFFLED) 

DF  =  8;r-cos(0J-^. 


(C-8) 


OCO  ^ 

x~M>{x,y)dxdy 


(C-9) 


DF  =  27:^cos(<P^). 

A‘ 

JJ  w(x,y)dxdy 

^  - = 

c  oo 

Jj  w\x,y)dxdy 


(C-10) 


(C-11) 


(C-14) 
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VOLUMETRIC  APERTURE 


\6k-  I  ,  7 


DF  =  — —-^m^m'^.cos'((l>^)  +  sin~((pJsm~(9J  +  rrr^ml  sin'((^,)cos^(0^) .  (C-15) 

A” 


m,  = 


:i:i 

“  ,  ^ 
x'w  {x,y,  z)dxdydz 

"OO 

V 

J^OQ  J--00 

[  w{x,y,z)dxdydz 

J— oo  y 

(C-16) 


DF  =  Ak^. 
A' 


(C-17) 


4  =  (w;4+^^4+-.^4r- 


(C-18) 


(\Wxsi^^y)d^d.y\ 

A^.  - —,  w„(x,y)  =  J  >v(x,y,z)Jz,  etc. 


\\wl.{x,y)dxdy 


(C-19) 


VOLUMETRIC  ARRAY 

\67r^ 


DF  =  — ^ ■yjm^m'^.cos^((pJ  +  mlm]  sm^(0Jsin^(0J  +  rtr^nF^  sin^(0Jcos^(0^)  .  (C-20) 

A“ 


V  —  2 


= 


V  '  J 


;  see  equation  (66). 


(C-21) 


DF  =  Ak^. 
A' 


(C-22) 


+v^,Nl,dld^  +  ulNl^ld]  . 


(C-23) 


r  Y 

_  V  _ J_ 


w 


^.{n,m)  =  ^w(n,m,k),  etc. 


(C-24) 
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OTHER  ARRAY  GEOMETRIES 


Directivity  approximations  for  many  additional  array  geometries  have  been  derived  in  the 
past/  DFs  for  some  other  common  types  are 


iTzdjX 


ring  array 


1  +  , 


K'd 


sm 


litd 


(C-25) 


where  d  is  the  diameter  of  the  continuous  ring  aperture.  For  djX  »  1 ,  the  directivity  factor 
approaches  -4  IndjX ,  which  increases  linearly  with  d.  Compare  this  behavior  with 

figures  10  and  1 1  for  a  semicircular  array. 


For  a  flat  circular  plate  of  diameter  d. 


DF 


{TtdIXf 


circular  plate 


1-- 


Tid  V  rrd 


sm 


iTtd  K 

~X 


(C-26) 


When  dlX»\,  this  equation  leads  to  the  well-known  approximation  for  a  piston  in  a  rigid  baffle, 
namely,  =  [tidl  A)  =  An  aJ  A' . 
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